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We present progress report of a CP-PACS calculation of quenched QCD spectrum with the Wilson quark 
action. Light hadron masses and meson decay constants are obtained at /3 —5.9, 6.1, and 6.25 on lattices with a 
physical extent of 3 fm, and for the range of quark mass corresponding to m^/mp ~ 0.75 — 0.4. Nucleon mass at 
each f3 appears to be a convex function of quark mass, and consequently the value at the physical quark mass is 
much smaller than previously thought. Hadron masses extrapolated to the continuum limit exhibits a significant 
deviation from experimental values: with K meson mass to fix strange quark mass, strange meson and baryon 
masses are systematically lower. Light quark masses determined from the axial Ward identity are shown to agree 
with those from perturbation theory in the continuum limit. Decay constants of mesons are also discussed. 



1. Introduction 

Precision determination of the light hadron 
spectrum is a fundamental problem of lattice 
QCD. Over the years numerous studies have been 
carried out on this subject, which have revealed 
a number of difficulties that have to be overcome 
to achieve this goal. In addition to the overrid- 
ing issue of how quenching distorts the spectrum, 
calculations aiming at precise results have to have 
good control over the systematic errors arising 
from a finite spatial lattice size, an extrapolation 
toward light quark masses, and an extrapolation 
to the continuum limit. There is also the issue 
of how reliably one can extract hadron masses 
from the propagators whose statistical errors gen- 
erally grow exponentially with the distance from 
the source point. 

Within quenched QCD high statistics calcula- 
tions on large lattices attempting to deal with 
these problems gradually developed since around 
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1990 exploiting the power of dedicated parallel 
computers In particular Weingarten and 

collaborators with the GF11 computer carried 
out a pioneering work on the continuum limit 
of the quenched spectrum with extensive simu- 
lations including finite size studies ||. In spite of 
these efforts convincing results on hadron masses 
in quenched QCD with a precision reliably better 
than 5% have not been obtained. In particular 
the question in what way the quenched spectrum 
deviates from the experiment has not been fully 
answered so far. 

In view of this situation the CP-PACS Collab- 
oration has decided to undertake a quenched sim- 
ulation of the light hadron spectrum as the first 
lattice QCD project of the CP-PACS computer 
which started operation in April 1996. The strat- 
egy we have chosen is to employ the simplest form 
of the action, i.e., the plaquette action for gluons 
and the Wilson action for quarks, and carry out, 
as much as the computing power of CP-PACS 
allows, a measurement of hadron masses down 
to small quark masses at small lattice spacings 
employing large lattices and high statistics. The 
simulation has been running since the summer of 
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1996 and are still continuing. In this article we 
report the present status of the work. 

2. Parameters of simulation 

In Table [j] we list the parameters chosen for 
our simulation, and the number of configurations 
for hadron mass measurements accumulated as of 
the time of the Workshop. 

We choose four values of the coupling constant 
(3 to cover the range of lattice spacing a -1 m 2 — 4 
GeV for controlling scaling violation effects and 
taking the continuum limit. In order to avoid fi- 
nite size effects, the spatial lattice size L is chosen 
so that the physical size equals approximately 3 
fm. The temporal lattice size is set equal to 7L/4; 
this choice is based on considerations on length of 
time intervals needed for a reliable extraction of 
masses with an exponentially smeared source we 
employ for hadron mass measurements. 

For each value of (3 we select five values of the 
hopping parameter such that the ratio m^/rrip 
takes values of about 0.75, 0.7, 0.6, 0.5, and 0.4. 
We abbreviate these hopping parameters as si, 
S2, Mi, U2 and 1*3. The first two values are for 
interpolation of results to strange quark, and the 
rest are for examining chiral extrapolation. Pre- 
vious spectrum studies with the Wilson quark ac- 
tion have been limited to the range m, v jm, p > 0.5. 
The value m^jm p ~ 0.4 represents our attempt 
toward lighter quark masses. Reducing the quark 
mass further is not easy: test runs we have carried 
out for m^/rrip « 0.3 at (3 = 5.9 show that fluctu- 
ations large and that it takes more computer time 
than that for the five hopping parameters down 
to m-n/mp 0.4. We calculate hadron masses for 
equal mass cases and also for the unequal mass 
cases of the type SiUj for mesons and SiSiUj and 
SiUjUj for baryons. 

Configurations are generated with the 5-hit 
Cabibbo-Marinari-Okawa heat-bath algorithm 
and the over-relaxation algorithm mixed in the 
ratio of 1:4. The combination is called a sweep, 
and we skip 200 to 2000 sweeps between hadron 
mass measurements as listed in Table [l]. 

Quark propagators are solved with the 
red/black-preconditioned minimal residual algo- 
rithm imposing the periodic boundary condition 



in all four directions. The stopping condition is 
chosen to ensure that truncation error in hadron 
propagators is at most 5 percent of our estimated 
final statistical error. For performance of our 
code on the CP-PACS see Refs. @§ 

Errors are estimated by a single elimination 
jackknife procedure. 

We note that results presented here are prelim- 
inary, as runs are still continuing. In particular 
the run at (3 — 6.47 has been just started. Al- 
though some figures include results at (3 — 6.47, 
values in the continuum limit are determined by 
extrapolating data at the first three values of (3. 

3. Extraction of hadron masses 

3.1. Smearing 

We employ an exponential smearing of quark 
source to enhance signals of ground state in 
hadron propagators. This choice is motivated 
by a recent measurement of the pion wave 
function defined by tp(r) — (0| J2 X ^( x )l5l{ x + 
r)k>/<0|£ B «(aOT69(a:)k) by the JLQCD Col- 
laboration, who found that ip(r) is well repro- 
duced by a single exponential function: ip{r) = 
Aexp(-5r)§. 

To choose the value of the slope B in our mea- 
surement, we parameterize their results in the 
formB(r7ip,m g ) = {x +xi-m p ) + (y +y 1 -(m p ) 2 )- 
m q . Estimates from available data for the value 
of m p appropriate for our simulation parameters 
are then substituted to find the value of B. The 
smearing radius 1 / B chosen in this way does not 
vary much over the points of our simulation, be- 
ing approximately given by 1/B « 0.33 fm. 
Various combinations of point (P) and smeared 
(S) sources are used to construct hadron opera- 
tors at the source. In the following we employ a 
notation such as PSS to specify the combination 
of quark sources. On the other hand, we always 
use point sinks for hadron operators at the sink. 

In Fig.|l| we show typical results of effective 
masses for various source combinations. We find 
that the effective mass reaches a plateau from 
above in almost all cases for all combinations of 
source. This suggests that the smearing radius 
of 1/B ?» 0.33 fm we choose is smaller than the 
actual spread of hadron wave functions. 
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Table 1 

Simulation parameters. Values of a -1 are determined from m p . 
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Figure 1. Typical results of effective masses 
obtained with various combinations of quark 
sources. 



we therefore decide to derive masses from the SS 
source for mesons and the SSS source for baryons. 

To extract hadron masses from propagators we 
employ a single hyperbolic cosine fit for mesons 
and a single exponential fit for baryons. The 
largest time slice of the fit i max is chosen by the 
requirement that the error of propagator does not 
exceed 5%. Changing the minimum time i m ; n of 
the fit, we compare results of correlated and un- 
correlated fits. We find that 1) x 2 /Nop for corre- 
lated fit decreases as t m [ n increases and becomes 
constant at some time slice called t x ^. 2) For 
£min > t x , mass results from correlated and un- 
correlated fits agree within at most 1.5 standard 
deviations. 3) Uncorrelated fits give results more 
stable as a function of t m i n . From these inves- 
tigations, we decide to adopt mass results from 
the uncorrelated fit at t m i n w t x , supplememted 
with a single elimination jackknife analysis to es- 
timate the error of masses. We consider that the 
magnitude of systematic error due to the choice of 
fitting ranges and fitting procedure is comparable 
to statistical error. 

In Table 0, we list typical statistical errors in 
mass results in units of percent, which are esti- 
mated by the single elimination jackknife method. 
In the contiuum limit results for the nucleon has 
the largest error of about 5% with present statis- 
tics. 



3.2. Fitting Procedure 

We observe in Fig. [j] that the onset of plateau 
is earliest if the smeared source is used for all 
quark and antiquark fields in hadron operators. 
Furthermore statistical errors of effective masses 
at a given time slice are the smallest for this com- 
bination. For results presented in this report, 



4. Extrapolation and interpolation to 
physical quark mass 

Our first interest is how hadron masses depend 
on the quark mass. We start with the pseudo 
scalar case, and plot in Fig.||the ratio m^/[(mi + 
7Ti2)/2] against (mi +m2)/2, where are quark 
masses defined by mi = Z m (l/Ki — 1/K C ). Con- 
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Table 2 

Errors in mass results in units of percent. 
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Figure 2. TO^/[(mi+TO 2 )/2] as a function of (mi + 
m 2 )/2. 

version to physical units is made with estimates 
of lattice spacing in terms of m p as described be- 
low. The critical hopping parameter K c is deter- 
mined by a linear extrapolation of m%, and Z m 
denotes the one-loop renormalization factor (see 
sec.^ for details). The ratio for 11 combinations 
of quark masses is constant within errors at each 
(3. Therefore we conclude that the pseudo scalar 
meson mass squared is linear in the average quark 
mass in the range covered by our data. In other 
words, we find no clear evidence for the existence 
of chiral perturbation theory higher order terms 
or quenched chiral logarithms^] in the results of 
pseudo scalar meson masses down to m q w 30 
MeV. 

Vector mesons and decuplet baryons also ex- 
hibit the property that their masses are linear in 
the average of quark masses with only a small 
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Figure 3. Pseudo-scalar meson masses squared 
and vector meson masses as functions of 1/K of 
light quark. 



mixture of higher order terms. This is illustrated 
in Fig|| and|I]in which hadron masses are plotted 
against 1/K of lighter quark, together with re- 
sults of a linear (pseudo scalar case) or quadratic 
(vector and decuplet case) fit. 

However, this property is badly broken for 
octet baryons as shown in Fig.[| Although S 
and 5 masses are linear in 1/K, nucleon and A 
masses are convex. The nucleon mass at the light- 
est quark mass is clearly off the linear fit (dashed 
line) as compared to a cubic fit shown by a solid 
line. 

In this connection we show our results for 
tojv/t7i p as a function of m^/rrip in Fig.| Dotted 
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Figure 4. Decuplet baryon masses against 1/K of 
light quark. 

and thick curves correspond to the cases in which 
the nucleon mass at (3 = 6.25 is fitted with a lin- 
ear or cubic functions of quark mass, respectively. 
Higher order terms are essential in order to make 
a reliable extrapolation of nucleon and A masses 
to the chiral limit, and have a significant impact 
on final results. 

Based on the considerations described above, 
we employ a linear fit for pseudo-scalar mesons, 
cubic fit for nucleon and quadratic fit for other 
hadrons in order to extrapolate or interpolate 
hadron mass results to physical quark masses at 
each j3. For vector mesons a linear fit reads to 
essentially the same result. We then determine 
values of K u ^ and a -1 taking the experimental 
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Figure 5. Octet baryon masses against 1/K of 
light quark. 

values of and m p . The hopping parameter 
K s for strange quark is estimated in two ways, 
either from assumed to be a pure ss state or 
from nix or my.. For the latter case we take 
meson mass data for unequal quark masses and 
first make an extrapolation of lighter quark to 
K u ,d with the heavier quark fixed at si and S2- 
We then linearly interpolate the result so that K 
or K* meson mass has the experimental value. 
For experimental values we use m T = 135 MeV, 
m p = 769 MeV, m = 1019 MeV, m K = 498 
MeV, and m K , = 896 MeV. 

5. Quenched hadron spectrum in the con- 
tinuum limit 



G 




240 



Figure 6. Edinburgh plot showing mjy/m p as a 
function of m v /m p . For explanation of lines see 
text. 

5.1. Scale parameter 

The scale parameter A is a fundamental quantity 
in QCD. We estimate the value in the MS scheme 
defined by 



A^=^(^(-/«))- 6l/26 "^ 1/26 °^ (7r/0) - 



where for the MS coupling we employ the 
tadpole-improved one-loop formula given by|lfj|, 
HI 



l/gj^(n/a) = tr(U P /3)/g 2 + 0.02461. 

We show the result in Fig.^|(a) as a function of 
m p a. Our results are consistent with those re- 
ported so far [§-||,[l| for the Wilson action. Our 
errors, however, are much reduced. In Fig.^(b) 
are compared our results with those for the 
Kogut-Susskind action|l2j and those from string 
tension E|. We find that all of the results agree 
in the continuum limit. 

5.2. Strange mesons 

In Fig.|| we show the continuum extrapolation 
for masses of mesons containing a strange quark. 
We find a linear behavior of the masses in lattice 
spacing to be well satisfied as shown by solid lines. 

As we noted in Sec. || strange quark mass can 
be determined in two ways. In Fig.||(a) we take 
ttik as input and predict the masses of cj> and K* 
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Figure 7. ^jjg vs. m p a (top figure) and a 1 
(bottom figure). Data for Wilson quarks are 
taken from GF11§, QCDPAX@, LANL 0, and 
APE||l[. Data for the KS action are from ref.|l^]. 
Data for the string tension are taken from ref. 1 13] . 

meson. The values extrapolated to the continuum 
disagree with experimental values plotted by open 
circles at a = 0. 

This should be contrasted with the results 
of the GF11 Collaboration^] plotted by open 
squares. They linearly extrapolate data measured 
on lattices with a spatial extent of about 2.3 fm 
(dashed lines). They then make a finite size cor- 
rection to the extrapolated value, and obtain a re- 
sult which agree well with experiment. The best 
fit to our data differ from those by GF11 taken 
on a 2.3 fm lattice, while their data on a 3.3 fm 
lattice at /3 = 5.7 are consistent with a linear fit 
of our results. 

We show an alternative procedure in Fig.||(b) 
where we plot results for rnx and ttik* when we 
take as input. In this case, the value of mx- 
in the continuum limit is consistent with experi- 
ment, while that for uik is much higher. These 
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Figure 8. Strange meson masses and their extrap- 
olations to the continuum limit. Open circles rep- 
resent experimental values. Data from GFll|| 
are plotted by open squares for comparison. 

results mean that there are no choice of strange 
quark mass which leads to an agreement of vector 
as well as pseudo-scalar meson masses. 

The J parameter [jl4j is a convenient quan- 
tity to summarize the situation. As shown in 
Fig.^|(a), the value of J in the continuum limit is 
smaller than the experimental value by more than 
2 standard deviations. In Fig.^(b), we compare 
our results with those obtained with the Kogut- 
Susskind quark actionjl^] and the clover quark 
action In the continuum limit, the two ac- 
tions give a result consistent with that of the Wil- 
son action. 

5.3. Baryon masses 

Our results for nucleon and A masses are shown 
in Fig.f[o|(a). For the nucleon our continuum ex- 
trapolated value is about 5% smaller than experi- 
ment. While the difference is only a one standard 
deviation effect, this contrasts with experience of 
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and refs. 15 Id]. We have calculated values of J 
using data of m n and m p in the cases that values 
of J are not reported in the literature. 



previous studies that the nucleon mass tends to 
come out higher than experiment. In fact our re- 
sults are significantly lower already at finite lat- 
tice spacings compared with previous results also 
shown in the figure. The discrepancy is mainly 
due to the bending of the nucleon mass toward 
small quark masses as discussed in Sec. ||. 

We emphasize that our results are consistent 
with those from other groups in the region of 
strange quark. If we make a linear chiral extrap- 
olation of the nucleon mass employing only three 
points of our data corresponding to m n /m p « 
0.75, 0.7, and 0.6, we obtain results in the chiral 
limit which agrees with those reported so far as 
illustrated in Fig.[To|(b). 

For A results from the present work and those 
of earlier studies overlap at finite lattice spacing. 
Nonetheless the continuum extrapolation leads to 
a value about 5% larger compared to experiment 
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Figure 10. (a) The same as Fig .|| for nucleon and 
A. (b) Nucleon masses obtained by linear fits to 
data at the largest three quark masses. 

in our case as compared to a value consistent with 
experiment reported by GF11. 

Figure |ll| shows results for strange baryons em- 
ploying m,K to determine the strange quark mass. 
We observe significant departures from the exper- 
iment except for £*. Our extrapolation is quite 
different from that of GF11 for fi and S* for which 
GF11 found an agreement with the experiment. 

5.4. Quenched spectrum of light hadrons 

We show all of our results for the quenched light 
hadron spectrum in Fig. [l| where we also plot the 
experimental spectrum (horizontal dashed bars) 
and the results of GF11 (open symbols) for com- 
parison. For strange quark mass the K meson 
mass is used as input both in our results and those 
of GF11. 

It is quite conspicuous that the masses of 
baryons containing strange quark are systemati- 
cally lower than experiment. This may be related 
to a similar deviation in the meson sector where 
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Figure 11. The same as FigJ8| for strange baryons. 



the masses of K * and <f> are smaller. 

It is interesting to note that agreement with ex- 
periment is much better if we examine the spec- 
trum within each flavor SU(3) multiplet. For ex- 
ample, if we take m p and to determine the 
lattice spacing and the strange quark mass, mass 
of K* is consistent with experiment, and a sim- 
ilar situation holds for E* and S* masses if we 
take the A and f2 masses instead for input. This 
property may contain some hint to consider the 
problem of how dynamical quark effects correct 



the discrepancy observed in Fig. 12. 

In the up-down quark sector our value for the 
nucleon mass is smaller than experiment by about 
5% and that for A larger by a similar amount. 
The errors are of a similar magnitude, however. 

Compared to the results of GF11, our errors are 
reduced by a factor two with our present statis- 
tics. The spectrum itself deviate from that of 
GF11 beyond estimated errors in a number of 
channels, most notable being cj), E* and fi. 
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Figure 12. Results for the quenched light hadron 
spectrum in the continuum limit from the present 
work as compared with experimental spectrum 
(horizontal dashed bars) and those of GF11 (open 
circles) . 

6. Light quark masses 

The values of light quark masses are important 
in a variety of phenomenological context |L7| . Per- 
turbatively quark mass is defined by 



2mP° rt - = Z m (l/A'! + 1/K 2 - 2/K c ) 



(1) 



with a renormalization factor Z m . A non- 
perturbative definition is also possible based on 
the axial vector Ward identity lia-pQ] 



(G|W|z) = 2m q (0\P\x). 

We take here the definition given by 



(2) 



2m 



PCAC = -m, ^ lim M^44, (3) 



where A4 is the local axial current and P is the 
pseudo-scalar density, and m^ on the righthand 
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Figure 13. Strange quark masses in the MS 
scheme at /1 = 2 GeV obtained from pertur- 
bation theory and from the axial Ward iden- 
tity. Rdiff(Ldiff) stands for right- (left-) differ- 
ence for an approximation to the time deriva- 
tive in the Ward identity. Values obtained with 
center-difference are numerically close to those 
with 771^, and are omitted in the figure. 

side is an approximation to the time derivative in 
the Ward identity. 

We compare results for the two definitions of 
quark mass. For the Ward identity method we 
test also e m " - 1, 1 - er m ^ and (e m - - e-' m -)/2 
instead of m, for replacing the time deriva- 
tive which correspond to right-, left- and center- 
difference, respectively. For the renormalization 
constants Z m j2]J], Za and Zp]2^|, we employ 
one-loop results improved by the tadpole pro- 
cedure with the coupling ay(l/a) [Io[| given by 
Z m = 8K C (1 + 0.01a v (l/a)), Z P = Z K (l - 
l.Q34a v (l/a)), and Z A = Z K (l-0.3Wa v (l/a)). 
Matching of lattice values of quark mass to those 
in the continuum is made at the scale 1/a em- 
ploying the MS scheme in the continuum. The 
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0.50 



Figure 14. The same as Fig.[13] for the average of 
up and down quark masses. 

result is then run by the two-loop renormalization 
group equation to 2 GeV. 

We plot results for the strange quark mass ob- 
tained with perturbative and three cases of non- 
perturbative PCAC definitions in Fig.|l3|. The <j) 
meson mass is used to determine K s in FigJT^(a) 
and the K meson mass is employed in (b). For 
finite lattice spacings the four definitions yield 
significantly different values. However, a linear 
extrapolation yields results in mutual agreement 
in the continuum limit. The same situation also 
holds for the average of up and down quark mass 



as shown in Fig. 14 



We note, however, that the value m s sa 140 
MeV for the strange quark mass in the continuum 
limit obtained with as input does not agree 
with m s ss 115 MeV found with thk as input. 
This is another sign of pathology in quenched 
QCD. 

In Fig.[l5], we compare our results for per- 
turbative quark masses with those from other 
groups ||,||-||]. Extrapolating our data to the 
continuum limit leads to m u .d ~ 4.2 MeV and 
m s « 136 MeV. They are larger than the val- 
ues estimated in Ref. |2^] [m u ^ ~ 3.4 MeV and 
m s » 96 MeV) based on a compilation of world 
data for the Wilson action. In the same figure we 
also plot results for the Kogut-Susskind action ob- 
tained by the JLQCD collaboration |l2] . It is not 
clear if Wilson and Kogut-Susskind results agree 
in the continuum limit. 
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Figure 15. Values of the average of up and down 
quark masses (top figure) and the strange quark 
mass (bottom figure). They are determined per- 
turbatively in the MS scheme at /z = 2 GeV. 
Data of hadron masses for Wilson quarks are 
taken from JLQCD§], GF11§, QCDPAX@ and 
LANL Q to calculate quark masses by the same 
method we employed for the CP-PACS data. Re- 
sults of quark masses for the KS action are taken 
from JLQCD|l2|. Solid squares at m p a = rep- 
resent best estimates in ref. 23 for the Wilson 
action in the continuum limit. 

7. Decay constants 

Determination of hadronic matrix elements of- 
ten suffers from uncertainties in renormalization 
factors. An exceptional case is the vector meson 
decay constant fv for which the renormalization 
constant for the local vector current V L can be 
precisely determined non-perturbatively with the 
use of the conserved current V c |Q . It has been 
noted by the QCDPAX Collaboration that this 
determination of the renormalization factor leads 
to a remarkable scaling of the decay constants Q]. 

In Fig.|l^ we plot /y obtained with the renor- 
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Figure 16. Results of f p and obtained with 
various renormalization constants. GF11 results 
|24[ are plotted for comparison. 

malization factor calculated non-perturbatively 
through Zy = 2KZy with Zy determined from 
(V C \V L )/(V L \V L ). Results obtained with naive 
perturbative Z factor to one loop given by Zy = 
2K (1 — 2.19 Qifatt.)p2| and with tadpole-improved 
perturbation theory Zy — (1 — 3K/4K C ) (1 — 
0.82 ay{\/a)) jl0];^5| are also shown for compar- 
ison. 

Making a linear extrapolation of the non- 
perturbative result in m p et, we obtain a value in 
the continuum limit which is completely consis- 
tent with experiment. Perturbative renormaliza- 
tion factors lead to slightly smaller values. It is 
worth noting that non-perturbative renormaliza- 
tion constant has a small value in the range of 
our data, given by Zy ss 0.54, 0.60, and 0.64 at 
[3 — 5.9, 6.1, and 6.25, respectively. A reliable ex- 
trapolation of the decay constants with one-loop 
evaluation of renormalization factors may be dif- 
ficult in such a situation. 

For pseudo-scalar meson decay constant fp, 



Figure 17. The same as Fig.|l6] for pseudo-scalar 
meson decay constants /„• and fx- 

the naive one-loop perturbative Z factor takes 
the form Z A = IK (1-1.68 a tott .) (22) and a 
tadpole- improved value is given by Za = (1 — 
3K/4K C ) (1 - 0.316 a v (l/a)) 0j2|. Both of 
these renormalization constants lead to a consis- 
tent value in the continuum limit as shown in 
Fig.[l7| Values of and Jk in the continuum 
limit are smaller than the experimental values by 
one to two standard deviations. 

A non-perturabtive estimate of the axial vector 
renormalization factor is available at /3 = 5.9 and 
6.l |26| , |27jl . Extrapolations of our data with these 
results for non-perturbative Za seems to give a 
consistent value with that with perturbative Za- 

8. Conclusions 

In this article we have presented a status re- 
port of our effort toward precision results of light 
hadron spectrum in quenched QCD with the stan- 
dard plaquette and Wilson quark actions. 

Our results have confirmed that strange quark 
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mass cannot be tuned in quenched QCD so that 
both pseudo scalar and vector spectra are in 
agreement with experiment. This observation 
is extended to the strange baryon sector in the 
present study, where our values for the baryon 
masses are systematically lower than experiment. 

We have also found that the nucleon mass de- 
creases faster than a linear behavior in quark 
mass toward the chiral limit. As a result the 
continuum value is smaller than the experiment, 
contrary to results of previous studies, albeit the 
difference is still a one standard deviation effect. 

Finally, we have shown that quark mass deter- 
mined from the axial Ward identity agrees with 
that from perturbation theory in the continuum 
limit, providing a handle for controlling the con- 
tinuum extrapolation. 

In parallel with these points of physics, an im- 
portant lesson we would like to draw from the 
present work is the importance of precision; with 
errors down to a percent level uncertainties in 
physics implications can be much removed. 

Encouraged by the results, we plan to further 
increase statistics of the three runs at (3 =5.9, 
6.1 and 6.25. In addition we shall soon have re- 
sults from the run at f3 = 6.47. We hope that 
results from the additional simulations help so- 
lidify our findings so far, and lead to a standard 
of the quenched hadron spectrum in QCD. 
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